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Abstract 

Many real systems possess accelerating statistics where the total number of edges 
grows faster than the network size. In this paper, we propose a simple weighted 
network model with accelerating growth. We derive analytical expressions for the 
evolutions and distributions for strength, degree, and weight, which are relevant to 
accelerating growth. We also find that accelerating growth determines the clustering 
coefficient of the networks. Interestingly, the distributions for strength, degree, and 
weight display a transition from scale-free to exponential form when the parameter 
with respect to accelerating growth increases from a small to large value. All the 
theoretical predictions are successfully contrasted with numerical simulations. 
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1 Introduction 



Standard interesting objects in the network science are relatively simple bi- 
nary (Boolean) networks where edges (links) are either present or absent, 
represented as binary states [T|2|3] . In other words, edges in Boolean networks 
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have equal "weights". However, the connections in many real networks are 
not homogeneous [I], which naturally calls for a typical measurement of the 
edge weight, such as the number of joint papers of two coauthors in scientific 
collaboration network [5,6,7,8], the number of flights or seats between any two 
cities in airline networks [9lfT0] . the bandwidth of a link in the Internet 
the reaction rate in metabolic network [12] , and so on. These real systems 
with diversity of edges can be better described in terms of weighted networks. 

Various weighted network models have been proposed to describe and explain 
the real- life systems [I] . Yook et al. took a first step in the direction of a model 
study for evolving weighted network (YJBT model) [13], where the topology 
and weight are driven by only the network connection based on preferential 
attachment (PA) rule [H]. The YJBT model overlooks the possible dynami- 
cal evolution and reinforcements of weights, which is a common property of 
real-life networks jl]. To better mimic the reality, Barrat, Barthelemy, and 
Vespignani presented a growing model (BBV) for weighted networks, where 
the evolutions of degree and weight are coupled [T5"fT6] . Enlightened by BBV's 
remarkable work, a variety of models and mechanisms for weighted networks 
have been proposed, including weight-driven model [T7], traffic-driven evolu- 
tion models [JSHH], fitness models [2D], local- world models [2"T|22[23] . deter- 
ministic models [2l|25ll26j . weight-dependent deactivation [27], spatial con- 
straints [281129] . 

Recent empirical study demonstrated that many real natural and social net- 
works exhibit the characteristic of "accelerating growth", which means that 
the total number of edges increases more quickly than linearly with the node 
(vertex) number. Generally, networks with this property are called "accel- 
erating networks". For instance, in both the World Wide Web [30J and the 
Internet [UJ the number of edges increases with time in a nonlinear fashion. 
In metabolic network [T2], the total number of links exceeds the total number 
of nodes by about an order of magnitude. Other familiar examples of acceler- 
ating networks include scientific collaboration network [5f6lT7ll8] . language net- 
work [31], citation network [32], etc. Inspired by this phenomenon, the research 
on accelerating networks has attracted an amount of attention [33| 34 35 3 6|37] . 

Actually, accelerating networks are far more common in the real world than has 
hitherto been appreciated [35J. Despite their widespread appearance, evolv- 
ing weighted networks with accelerating growth have received less attention. 
This important factor of accelerating growth is neglected in most of the con- 
Then questions arise naturally: Can an accelerating growth model for weighted 
networks be established? How much effect does the accelerating growth have 
on the evolution of weighted networks? 

In this paper, we present an accelerating weighted network model to under- 
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stand how the accelerating growth phenomenon affects the dynamical evolu- 
tion of weighted networks. We study both analytically and numerically the 
network characteristics, including the evolution and distributions of the de- 
gree, weight, and strength, as well as the clustering coefficient. We show that 
obtained properties depend on the accelerating parameter. 



2 Definitions 

We give a brief introduction to the definitions of tools for statistical charac- 
terization of weighted networks. 

Mathematically, for a weighted network, its topological as well as weighted 
properties can be completely described by a generalized adjacency matrix W, 
whose element Wij specifies the weight of the edge between node i and j. 
Wij = represents that node % and j are disconnected. In the following, we 
focus on the cases of undirected graphs, which have symmetric nonnegative 
weights = Wji > 0. Moreover, we assure that Wa = 0. 

The standard topological characterization of binary networks is also applied for 
weighted networks, which is obtained by the analysis of the distribution P(k) 
that represents the probability of a random selected node to have degree k. 
In a weighted network, a natural generalization of degree is the node strength 
defined as Sj = X^er(i) w ij-> where the sum runs over the set F(i) of neighbors 
of the node i. Statistical properties of weighted networks can be characterized 
by the distributions of strength P(s) and weight P(w), which denote the 
probability of a node to have strength s and of an edge to have weight w. 



3 The model for accelerating weighted networks 

The model proposed here begins from an initial configuration of Nq nodes con- 
nected by edges with assigned weight wo — 1. At each time step, the network 
evolves under the following two coupled mechanisms: topological growth and 
weights' dynamics. 

(i) Topological growth. A new node n enters the network. If the new node is 
born at time t, we assume that the edge number of this new node is a power 
law function t 6 (0 < 9 < 1) depending on time t, where 9 is the acceleration 
parameter. (Since multiple edges are forbidden, the total number of edges is 
smaller than t 2 /2, thus 9 cannot be greater than 1. On the other hand, when 
6 < 0, the new node may carry no edge, which is not consistent with most 
real networks. Thus, one may reasonably assume < 9 < 1.) These t e new 
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edges have initial weight wq = 1 and are randomly attached to a previously 
existing node i according to the preferential probability 



n 



-I 



J2j Sj 



(1) 



where the sum runs over all existing nodes. 

(ii) Weights' dynamics. The creation of each of the t e edges will introduce 
variations of the existing weights across the network. For the sake of the sim- 
plicity, we only consider the local rearrangements of weights of those edges 
connecting i and its neighbors j e according to the simple rule 

w i: j -> Wij + 5— . (2) 



Here we have assumed the addition of each new edge induces a total increment 
5 (5=const) of weights. The rule described by Eq. (T5|) yields a global increase 
of Wo + S for the strength of node i, which will therefore become even more 
attractive to future nodes. 

After updating the weights, the growing process is iterated by introducing 
another new node, i.e. returning to step (i) until the network reaches the 
desired size. Since the network size is incremented by one with each time step, 
we use the step value t to represent the node created at this step. At time 
t, the network has N = t + N nodes and / t e dt = t 1+9 /(l + 6) edges in the 
continuum limit. 

Note that many real-life networks exhibit such an evolving mechanism as de- 
scribed in our model. Mechanism (i) is a plausible one that appears in many 
real systems such as the Internet and scientific collaboration networks; it cor- 
responds to the fact that new nodes try to connect a preexisting node with a 
probability proportional to the strength of the old node, as provided by Eq. JT]). 
An important aspect of mechanism (i) is that the number of edges carried by 
the new node is not a constant, but controlled by a parameter 6 (0 < 6 < 1), 
which we call accelerating exponent. This has been confirmed by a variety 
of empirical observations. For example, for arXiv citation graph, autonomous 
system graph of the Internet, and the email networks, their accelerating ex- 
ponents have been found to be 0.56, 0.18, and 0.12, respectively [38J. On the 
other hand, mechanism (ii) describes the weight dynamics induced by a new 
edge onto the old ones, which can happen in scientific collaboration networks, 
airline networks, and so on. See Refs. [9|15II16] for detailed explanation. 




The model is governed by two parameters 5 and 6, according to which there 
are some limiting cases of the model. When 5 = and 6 = 0, it is reduced 
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to the BA model [H]. For 9 = 0, it coincides with a special case of the BBV 
model [ToTTE] . Varying 5 and 9 allows one to study the crossover between the 
two limiting models, which have qualitatively different properties from the 
two limits. We will show that both of the parameters 5 and 9 have significant 
effects on the network evolution, here we focus on the latter that has never 
been studied before, while the former has been discussed in Ref. |16j . 



4 Evolution and distributions of strength, degree and weight 



Our growing model can be studied analytically through the time evolution 
of the average value of Si(t) and kt(t) of the ith node at time t by using the 
"mean-field" method. According to the evolution rules, the addition of each 
new edge results in the increase of total network strength by an amount equal 
to 2 + 25. Thus, after t steps of evolution, the total strength of the network is 
obtained to be 



$>(f) =£A2 + 25) 

i=0 

t 

i\2 + 25)di 

] -t e+1 . (3) 



o 



2(1 + 



9 + 1 



When a new node n enters the network, an existing node i can be affected in 
two ways: (1) The new node is connected to i with probability given by Eq. 
([TJ), thus the degree and strength of i are increased by 1 and 1 + 5, respectively. 
(2) The new node is connected to one of i's neighbors j G in this case the 
degree of % remains unmodified, while tUy is increased according to Eq. (T5]), 
and thus S{ is increased by 5wij/sj. We assume that variables Si and ki are 
continuous. Then, at each time step, the strength s» and degree ki of a node 
evolve as 



dsj 
dt 



and 



dki 
~dt 



Sj ^j 



25 + 1 
25 + 2 [ 



Si(t) 



9 Si (t) 



> Sj(t) 

J2jSj{t) 2 + 25 t 
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respectively. 



With the initial conditions ki(t = i) = Si(t = i) = i e , we can integrate above 
equations to obtain 

R (25+l)(l+9) 



and 



ki(t) 



(!) 



Si{t) + 25t d 
25+1 

(25+l)(l+g) 

25+2 + 25t e 

—■ (7) 



25+1 

For large t, it can be always guaranteed that Si(t) is larger than 25t e because 
one can easily see that ~^+2 + ^ > ^ gi ven $ G [0, 1] and 5 > 0. Thus, Sj(i) 
and fcj(t) are proportional when t is large. From the obtained expressions, 
we find that dynamical exponent (3 = ^ 2<5 ^^ +e - > depends on not only S, but 
also 9, which is significantly different from the BBV model and BA network. 
Figure [1] shows the behavior of the nodes' degree and strength versus time for 
different 9, which recovers the results predicted analytically. 

We assume that the nodes are added to the systems at equal time intervals, 
then the probability that a node has strength Si(t) smaller than s, P(si(t) < s), 
can be written as 



. y \ . „/ 25+2 (25+l)(l+fl) \ 

P{Si(t) < S) — P U > S 25+i-et 24+i-e 1 



25+2 (25+l)(l+9) 
1— P [i<S 2S+l-9f 25+1-0 



1 25+2 (25+l)(l+9) 

= 1 S 25+1-flt 25+1-6 , (8) 

The probability distribution of strength P(s) can be calculated through solving 
the partial differentiation of P(si(t) < s) on s, and the final result of strength 
distribution at time t exhibits the following behavior: 



where f(t) = _L_ t( 2<5 + 1 )( 1 + 6 ')/( 2<5 + 1 e ) . Therefore, the strength follows a power 
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Fig. 1. (Color online) Time-evolution for the degree and strength of nodes, added 
to the systems started at t = 1 in the case of 5 = 2. 

law distribution with an exponent 



Is 



45 + 3-6* 



25 



9 ' 



(10) 



Since there is an approximatively linear relationship between the degree and 
strength for the same node, the degree distribution P(k) has also a scale-free 
form, P(k) ~ k~ lk with the exponent 7& identical to 7 S . 

In order to confirm the validity of the obtained analytical predictions, we per- 
formed extensive numerical simulations of the networks. To reduce the effect 
of fluctuation on simulation results, the simulation results are average over ten 
network realizations. Figure [2] gives the accumulative distributions of nodes' 
degree and strength. Numerical results are consistent with the theoretical ones. 

Now we investigate the evolution of the weights with time, which can also 
be computed analytically using mean field approximation employed for the 
research of s*(t) and ki(t). During the process of network growth, Wij can only 
increase by the addition of a new node connected to either node i or j, and 
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Fig. 2. (Color online) The cumulative distributions of nodes' degree and strength 
at various 9 values for 5 = 5. The network size is N = 5000, and the straight lines 
are the theoretical results of ^k^cum = 1s,cum = 

(4J + 3 - 0)/(2tf + 1 - 0) - 1. 
the evolution of satisfies the following equation 



dwij 



J2k s k Si 



+ t< 



J2k S k Sj 



6{l + 6) w i:i 
1 + 5 t 



The edge is built only when both node i and j have been created, there- 
fore the birth time of edge is = max(i,j). Considering the initial 
condition Wij = 1, one can integrate the Eq. (1111) to obtain 

6(1+6) 
{ t \ 6+1 

^ = [ J ' (12) 



which implies that the edge weight w displays a power-law distribution, P(w) ~ 
w~ lw , with the exponent j w = ^n^jjf- ■ To check the validity of the analyti- 
cal predictions for the evolution and distribution of the edge weight, we have 
performed numerical simulations of the present model, which are plotted in 
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Fig. 3. (Color online) (Upper panel) Log-log plot of time-evolution for the weight of 
edges, which are born at step t = 1 in the case of 5 = 2. (Lower panel) The cumula- 
tive distribution of edge weight for network size N = 5000 at different 6 and 5. The 
dashed straight lines are the theoretical results of 7 tojCum = (1+25+59) /[(l+6)5] — 1. 



Fig. [31 The simulations are in agreement with the analytical calculations. 



So far, we have shown that the considered accelerating weighted network has 
power-law distributions of strength, degree, and weight. All the obtained ex- 
ponents 7fc = 7s = H+fzf j an d lw = 1 (i+g)g 6> vary from 2 to oo, depending 
on the network parameters 5 and 9. It should be noted that in the BBV 
model, where 9 equals zero, exponents 7^ = 7 S are between 2 and 3. Thus, 
accelerating growth has an important effect on the evolution of the network: 
when 9 increases from to 1, distributions of strength, degree, and weight 
exhibit a transition from scale-free (small exponent) to exponential (large ex- 
ponent) forms. Moreover, by tuning the values of parameter 9, the network 
may have different forms of degree (strength) distribution and weight distri- 
bution. For example, when 5 — * 0, the exponents 7^ (7 S ) and 7^ are reduced 
to 7^ = 7 S = 1 + yzg and 7^ = 00, respectively. In this case, the weight 
distribution is always exponential, while the degree (strength) distribution 
follows a power-law form with 7^ (7,,) increasing from 3 to 00 when parameter 
9 increases from to 1. 
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Fig. 4. (Color online) Average clustering coefficient C vs parameter 9 and the net- 
work size N . All data are from the average of ten independent simulations for the 
same network. 

5 Clustering coefficient 

As studied in the previous section, accelerating growth significantly affects the 
network properties, such as distributions of strength, degree, and weight. In 
this section, we will show that the parameter 6 concerning accelerating growth 
also controls the clustering coefficient of the networks. 

By definition, the clustering coefficient [39] of a given node is the ratio of the 
total number of edges that actually exist between all its k nearest neighbors 
and the potential number of edges k(k — l)/2 between them. The clustering 
coefficient of the whole network is obtained through averaging over all its 
nodes. 

We have performed numerical simulations of the networks to study the in- 
fluence of the acceleration parameter 6 on clustering coefficient, which is 
presented in Fig. HI Numerical simulation results show that for arbitrary 6, 
the average clustering coefficient C(N) decreases as the size increases, i.e. 
C(N) ~ N~ v , as reported in Fig. HI This power-law behavior is similar to 
that of Barabasi- Albert (BA) model pQ, and is in contrast to the observa- 
tion of some real networks such as the Internet and metabolic networks whose 
clustering coefficient is independent of their size [ID]. However, to the best of 
our knowledge, whether there is real networks exhibiting similar phenomenon 
of clustering coefficient as our model is still unknown, which deserves further 
study in future. 

In fact, there are two known limits for the average clustering coefficient C of 
the whole network: In the case of 6 = 0, the network is a tree, and hence it 
has no triangle, therefore C = 0; when 6 = 1, the network corresponds to a 
complete graph (i.e. A^-clique) with clustering coefficient C = 1. Figure H] also 
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shows the dependence of C on 6. One can easily see that C is an increasing 
function of 9 as expected. The increase is not very sharp for small 6, but we 
can expect an obvious increase of C for 9 of very large values. 



6 Conclusion 



In conclusion, we have proposed a growing model for weighted networks in 
which the number of edges added with each new node is an increasing function 
(power law function) of network size. We have demonstrated that the accel- 
erating growth is an important factor that establishes the network structure. 
Using mean field network theory, we have computed analytical expressions for 
the evolution and distributions of strength, degree and weight. The obtained 
results show that these distributions are subject to a transition from a power- 
law to exponential shape, when the acceleration parameter is tuned from a 
small to large value. All mean field approximations have been confirmed by nu- 
merical simulations. Our model may provide valuable insight into the real-life 
networks. 

Although accelerating growth exists in many real-life networks, it should be 
pointed out that the nonlinear growth fashion of the edges compared with 
nodes in real systems is more intricate and flexible. We use here the most 
generic case, i.e. the number of total edges is a power-law function of the total 
node number. In future, it would be worth studying in detail other manners 
of nonlinear growth in different real-life networks as well as their impacts on 
network properties and dynamics. 



Acknowledgments 

The authors are grateful to the anonymous referees for their valuable com- 
ments and suggestions. This research was supported by the National Basic 
Research Program of China under grant No. 2007CB3 10806, the National 
Natural Science Foundation of China under Grant Nos. 60496327, 60573183, 
60773123, and 60704044, the Shanghai Natural Science Foundation under 
Grant No. 06ZR14013, the Postdoctoral Science Foundation of China un- 
der Grant No. 20060400162, Shanghai Leading Academic Discipline Project 
No. B114, the Program for New Century Excellent Talents in University of 
China (NCET-06-0376), and the Huawei Foundation of Science and Technol- 
ogy (YJCB2007031IN). 



11 



References 

[I] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 (2002). 

[2] S. N. Dorogvtsev and J.F.F. Mendes, Adv. Phys. 51, 1079 (2002). 
[3] M. E. J. Newman, SIAM Review 45, 167 (2003). 

[4] S. Boccaletti, V. Latora, Y. Moreno, M. Chavezf, and D.-U. Hwanga, Phy. Rep. 
424, 175 (2006). 

[5] M. E. J. Newman, Proc. Natl. Acad. Sci. U.S.A. 98, 404 (2001). 
[6] M. E. J. Newman, Phys. Rev. E 64, 016132 (2001). 

[7] A.-L. Barabasi, H. Jeong, Z. Neda. E. Ravasz, A. Schubert, and T. Vicsek, 
Physica A 311, 590 (2002). 

[8] M. Li, J. Wu, D. Wang, T. Zhou, Z. Di, Y. Fan, Physica A 375, 355 (2007). 

[9] A. Barrat, M. Barthelemy, R. Pastor-Satorras, and A. Vespignani, Proc. Natl. 
Acad. Sci. U.S.A. 101, 3747 (2004). 

[10] W. Li, and X. Cai, Phys. Rev. E 69, 046106 (2004). 

[II] M. Faloutsos, P. Faloutsos and C. Faloutsos, Comput. Commun. Rev. 29, 251 
(1999). 

[12] H. Jeong, B. Tombor, R. Albert, Z. N. Oltvai and A.-L. Barabasi, Nature 407, 
651 (2000). 

[13] S. H. Yook, H. Jeong, A.-L. Barabasi, Y. Tu, Phys. Rev. Lett. 86, 5835 (2001). 
[14] A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

[15] A. Barrat, M. Barthelemy, and A. Vespignani, Phys. Rev. Lett. 92, 228701 
(2004). 

[16] A. Barrat, M. Barthelemy, and A. Vespignani, Phys. Rev. E 70, 066149 (2004). 

[17] T. Antal and P. L. Krapivsky, Phys. Rev. E 71 026103 (2005). 

[18] W.-X. Wang, B.-H. Wang, B. Hu, G. Yan, and Q. Ou, Phys. Rev. Lett. 94, 
188702 (2005). 

[19] Y.-B. Xie, W.-X. Wang, and B.-H. Wang, Phys. Rev. E 75, 026111 (2007). 

[20] G. Bianconi, Europhys. Lett. 71, 1029 (2005). 

[21] Z. Pan, X. Li, and X. Wang, Phys. Rev. E 73, 056109 (2006). 

[22] B. Wang, H. W. Tang, Z. Z. Zhang, and Z. L. Xiu, Int. J. Mod. Phys. B 19, 
3951 (2005). 



12 



[23] Z.Z. Zhang, L. L. Rong, B. Wang, S. G. Zhou, and J. H. Guan, Physica A 380, 
639 (2007). 

[24] S. N. Dorogvtsev and J. F. F. Mendes, AIP Conf. Proc. 776, 29 (2005). 

[25] Z.Z. Zhang, S.G. Zhou, L.J. Fang, J.H. Guan, Y.C. Zhang, EPL 79, 38007 
(2007) 

[26] Z.Z. Zhang, S. G. Zhou, L. C. Chen, J. H. Guan, L. J. Fang, and Y. C. Zhang, 
Eur. Phys. J. B 59, 99 (2007). 

[27] Z.-X. Wu, X.-J. Xu, and Y.-H. Wang, Phys. Rev. E 71, 066124 (2005). 

[28] A. Barrat, M. Barthelemy, and A. Vespignani, J. Stat. Mech: Theory Exp. 
P05003 (2005). 

[29] G. Mukherjee and S. S. Manna, Phys. Rev. E 74, 036111 (2006). 

[30] A. Broder, R. Kumar, F. Maghoul, P. Raghavan, S. Rajagopalan, R. Stata, A. 
Tomkins, and J. Wiener, Comput. Netw. 33, 309 (2000). 

[31] S. N. Dorogovtsev and J. F. F. Mendes, Proc. R. Soc. London, Ser. B 268, 2261 
(2001) 

[32] S. Redner, Eur. Phys. J. B 4, 131 (1998). 

[33] S. N. Dorogovtsev and J. F. F. Mendes, Phys. Rev. E 63, 025101(R) (2001). 

[34] P. Sen, Phys. Rev. E 69, 046107 (2004). 

[35] J. S. Mattick and M. J. Gagen, Science 307, 856 (2005). 

[36] M. J. Gagen and J. S. Mattick, Phys. Rev. E 72, 016123 (2005). 

[37] X. Yu, Z. Li, D. Zhang, F. Liang, X. Wang, and X. Wu, J. Phys. A 39, 14343 
(2006). 

[38] J. Leskovec, J. Kleinberg, C. Faloutsos, ACM Transactions on Knowledge 
Discovery from Data, 1, 1 (2007). 

[39] D. J. Watts and H. Strogatz, Nature (London) 393, 440 (1998). 

[40] E. Ravasz, A.-L. Barabasi, Phys. Rev. E 67, 026112 (2003). 



13 



